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The International Gravitational Event Collaboration, IGEC, is a coordinated effort by 
research groups operating gravitational wave detectors working towards the detection of 
millisecond bursts of gravitational waves. Here we report on the current IGEC resonant 
bar observatory, its data analysis procedures, the main properties of the first exchanged 
data set. Even though the available data set is not complete, in the years 1997 and 1998 
up to four detectors were operating simultaneously. Preliminary results arc mentioned. 

1. The International Gravitational Event Collaboration 

One of the most relevant scientific objectives for resonant detectors of gravitational 
waves (gw) is the search for short gw bursts, emitted during the gravitational col- 
lapse of stars or the final evolution of coalescing binaries .□ To ensure the confidence 
of a detection, it is necessary to compare the observations made by multiple detec- 
tors with uncorrelated noise. This has already been done in the past years among 
pairs of crYJDgenic detectors with common time periods of observation of about a 
semester .B'ta A few days of observation have been reported also for three simulta- 
neously operating detectors.H In these attempts, the gw search consisted of a time 
coincidence analysis among the candidate signals reported by the different detectors 
and no statistically significant excess of coincidences was found. 

An increase of the number of cryogenic resonant detectors in simultaneous opera- 
tion in recent years has greatly improved the chance of making a confident detection 
of gw bursts. In fact, the International Gravitational Event Collaboration, IGEC, 
currently consists of the research groups operating the five cryogenic bar detectors 
ALLEGRO,! AURIGA,! EXPLORER^! NAUTILUSQ and NIOBE.! The IGEC was 
established in July 1997 with an agreement^ for setting up a common search for 
gravitational wave bursts of duration of the order of 1ms. This agreement sets the 
guidelines for the data exchange procedure among the participating groups and the 
IGEC scientific policy, whose most relevant aspects are: 

• each group has responsibility to make available to IGEC its list of candidate 
gravitational wave events, 

• a unanimous agreement of the member groups is required to make public the 
results based on the IGEC data exchange, 

• IGEC is open to new data taking research groups. 

In 1999 the first IGEC analysis of the 1997-1998 data was performed and some initial 
results will be presented in the following sections. At the time of this analysis, not 
all the 1997-1998 data had been exchanged. Despite the incomplete data set, the 
simultaneous operation of four gravitational wave detectors was achieved for the 
first time. 
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1.1. Data exchange protocol 

The IGEC data exchange procedure is aimed at searching for coincident excitations 
at different detectors. For each detector, a list of candidate events, each describing a 
5-like gravitational wave excitation of the detector, is provided by the corresponding 
research group. The IGEC protocol requires that in each list the candidate event 
rate be at most of the order of 100 /day, to limit the expected rate of accidental co- 
incidences. Currently, the research groups have been exchanging event lists relative 
to the past three years, but a future goal is to implement an automatic exchange 
day by day. 

The event lists^re then made available to the IGEC collaboration as files under a 
common protocol,™ open for future extensions. Under this protocol, it is mandatory 
for each detector to provide the minimum set of information needed to describe 
a (5-like gw excitation of the detector for each event; namely its Universal Time 
of arrival, the Fourier component H Q of its amplitude in Hz -1 and the detector 
noise level at that time. Another mandatory requirement is the declaration of the 
effective observation time of each detector, so that the IGEC observation time can 
be calculated. Optional information, such as the time of threshold crossing of the 
detector output, the duration of the event and its statistical compliance to a (5-like 
gw excitation, can also be exchanged. 

1.2. IGEC gravitational wave observatory 

The relevant parameters of the five resonant bar detectors of the IGEC observatory 
in the years 1997-1998 are summarised in Table 1. The detectors are sensitive 
to gw signals in a typical bandwidth of the order of 1Hz around each one of the 
two resonances of the detector, which are close to 700Hz for NIOBE and close 
to 900Hz for all other detectors. The relationship between the Fourier amplitude 
H , averaged on the two resonant frequencies of the detector, and the energy E s 
deposited by the g.w. burst on the bar, is given by: 

H = - T 1 2 y/Es/Mbar (1) 

where L\, ar is the bar length, Mt, ar its mass, vq the mean of the detector resonance 
frequencies. 

The typical thresholds used for selecting the events in the 1997-1998 burst search 
have been in the range Hq — 2 — 6 x 10~ 21 Hz^ 1 . The corresponding strain am- 
plitude of the gw can be computed assuming a model for the burst shape: for 
the conventional ~ 1ms burst, the Fourier component Hq should be multiplied by 
10 3 Hz to get the maximum strain amplitude h. 

To maximise the chances of a coincidence detection, the bars have been oriented 
to be approximately parallel to one another. Neglecting the polarisation effects, 
the gw amplitude at the detector is Hq = H gw sin 2 9(t), where H gw is the incident 
gw amplitude and 9{t) is the angle between the bar axis and the direction of the 
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Table 1. Main characteristics of the five resonant bar detectors in 1997-1998. The reported 
misalignment is the angle between the bar axis and an optimal direction whose overall misalignment 
from the detectors is minimal. Q± is the typical quality factor of the resonances 



ANTENNA 


ALLEGRO 


AURIGA 


EXPLORER 


NAUTILUS 


NIOBE 


Material 


AL5056 


AL5056 


AL5056 


AL5056 


Nb 


Mass [kg] 


2296 


2230 


2270 


2260 


1500 


Length [m] 


3.0 


2.9 


3.0 


3.0 


2.75 


Mode - [Hz] 


895 


912 


905 


908 


694 


Mode + [Hz] 


920 


930 


921 


924 


713 


Q± [io 6 ] 


2 


3 


1.5 


0.5 


20 


Temp. [K] 


4.2 


0.2 


2.6 


0.1 


5.0 


Longitude 


268°50'£ 


11°56'54"E 


6°12'£! 


12°40'21"E 


115°49'E 


Latitude 


30° 27' iV 


45°21'12"2V 


46°27'A r 


41° 49' 26" AT 


31°56'5 


Azimuth 


40° W 


44° E 


39° E 


44° E 


0° 


Misalignment [deg] 


9 


4 


2 


3 


29 



source. In this configuration of the observatory, the relative misalignments reported 
in Tab. 1 for ALLEGRO, AURIGA, EXPLORER and NAUTILUS disperse their 
sin 2 9(t) responses by at most a few %, while for NIOBE the dispersion is up to 
a few tenths. Figure 1 shows as an example the resulting amplitude efficiency for 
the observation of the Galactic Center as the Earth rotates. Since their values of 
sin 2 6(t) are simultaneously above 0.7 for about 60% of the time, we point out that 
this configuration of the observatory ensures a rather good and coherent coverage 
of the central galactic mass during time. 




Time (h) 



Fig. 1. Amplitude response sin 2 9(t) of the five IGEC detectors versus Universal Time for signals 
incoming from the Galactic Center during one day. 



2. The Data Set 

Each IGEC group independently developed a data acquisition and an optimum 
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filtering procedure for a 5-like gw excitation. These procedures differ greatly in 
the methods^in particular some of them filter for the energy released in the bar 
(ALLEGRO ,□ NIOBEEI) while the others filter for the amplitude and phase of the 
strain excitation of the bar (AURIGaJLj EXPLORER and NAUTILUS0). The 
output correlation time of the filters ranges from a few tenths up to a few seconds. 

A search for maxima on the filter outputs is then used to identify the time and 
amplitude of the candidate events, which are exchanged only if their amplitude ex- 
ceeds a selected threshold relative to the noise level. For the current detectors, these 
thresholds span in the range of signal-to-noise ratio SNR ~ 3 — 5 in amplitude. 
Of the events overcoming the threshold, some are rejected as spurious by different 
methods. AURIGA implements a \ 2 test on the conrrjiance of the single detected 
excitations with the expected template of a gw burst S3 In fact, the filtering proce- 
dure implemented is equivalent to a maximum likelihood fit of a signal model to the 
data and the goodness of the fit is statistically tested for each candidate event; those 
events not passing the test are rejected. All the other detectors implement spurious 
signals rejection using sensors of ambient disturbances. In addition, EXPLORER 
and NAUTILUS reject the events when the local noise is above a certain threshold. 

The effective observation times for each detector has been defined by vetoing, a 
priori, time periods of detector maintenance or malfunctions and times when the 
detector was excited by the local environment as determined by the local experimen- 
talists. After filtering, AURIGA implements a second level of vetoes a posteriori to 
reject further periods of unsatisfactory performance due to a lack of self consistency 
of its data analysis, that is when its noise fails to be compliant with the modeled 
one used to build its filtering procedures 



Table 2. Left: total observation time T ^ s and rate of events R ev t for each detector. Right: net 
common observation time Tjv when at least N^etectors were simultaneously operating. 



Detector 


T ohs (day) 


Revt (events / day) 


■*V detector s 


T N (day) 


ALLEGRO 


405.7 


112.9 


l 


625.0 


AURIGA 


153.0 


175.3 


2 


260.4 


EXPLORER 


137.5 


150.7 


3 


89.7 


NAUTILUS 


108.5 


80.8 


4 


15.5 


NIOBE 


185.9 


14.0 


5 






The amplitude distributions of all the 1997-1998 exchanged events for each de- 
tector are shown in Fig. 2. The effective observation time of the exchanged data up 
to now is summarised in Table 2, together with the mean rate of exchanged events. 
The net observation time with at least four, three and two detectors simultane- 
ously operating has been respectively 15.5, 90 and 260 days. We expect that the 
three-way observation time will increase significantly as the exchanged data set will 
become complete. The ALLEGRO detector has been showing the best duty cycle, 
close to 100% on the exchanged data period, as well as the most stationary noise 
performance with respect to the other detectors. 
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Fig. 2. Density functions of the amplitude of all the exchanged events in 1997-1998 for each detec- 
tor. The densities are normalised to unit area and are estimated from the amplitude histograms 
of the events. 



3. Analysis of Time Coincidences 

A search for two, three and four-fold coincidences was carried out on the exchanged 
data. In the analysis reported here a M-fold coincidence is observed if the estimated 
arrival times U at the M detectors are all — tj\ < Is. This figure has been chosen 
as a compromise between the demands for a small accidental background and for 
a low false dismissal probability.-'' In fact, the measured uncertainties t w on the 
estimated arrival times of a burst at each detector are quite similar and are selected 
to be ±0.5s. This corresponds to a maximum false dismissal of a few % for the 
exchanged events, even at low SNR. The maximum separation among coincident 
events is therefore set to 2t w = Is. 

A preliminary search for three and four-fold coincidences shows none. A detailed 
analysis is in progress. The two-fold coincidences found for each pair of detectors 
are shown in Table 3: in all the cases they correspond to the estimated accidental 
background. 

Two methods for estimating the rate of accidental coincidences have been applied 
for the pairs of detectors: i) performing several time shifts of events times of one 
detector with respect to the other and then looking for coincidences! (the resulting 
accidental coincidences are reported as < n a > in Tab. 3); ii) assuming stationary 
Poisson distributions of event times and using the mean measured rates of events 
for each detector during the common observation time. The latter method predicts 

t a new approach for multiple time coincidence analysis has been proposed by the Rome group.Q 
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Table 3. Preliminary results of two way coincidence analysis for the 1997-1998 data. The ab- 
breviations, AL, AU, EX, NA and NI stand for ALLEGRO, AURIGA, EXPLORER NAUTILUS 
and NIOBE respectively. n c is the found number of coincidences, < n a > and < n a >theory the 
expected accidental ones respectively from the time shift method and Eq.^, P(n > n c \ < n a >) 
the calculated probability for the coincidences to be > n c , and T (, s the common observation time 
of the pairs of detectors. 



Detectors 


n c 


< n a > 


Tl>a ^theory 


P(n > n c \ < n a >) 


^obs 


AL-AU 


42 


46.6 


45.1 


0.77 


103.8 


AL-EX 


27 


31.2 


30.9 


0.80 


100.7 


AL-NA 


17 


21.6 


21.0 


0.87 


98.9 


AL-NI 


1 


0.9 


1.0 


0.61 


27.1 


AU-EX 


14 


20.3 


19.2 


0.94 


44.1 


AU-NA 


4 


4.2 


4.2 


0.60 


18.3 


AU-NI 


1 


2.3 


2.1 


0.90 


37.0 


EX-NA 


5 


7.0 


5.7 


0.83 


37.5 


EX-NI 


1 


1.1 


1.0 


0.65 


18.9 



a number of accidental M-fold coincidences given byEi 

< n a >theory= M " _ 1 71*, (2) 

1 obs i=l 

where M is the number of detectors, T b s their common observation time, t w the 
time window describing the timing accuracy of each detector, n, the number of 
exchanged candidate events of the i th detector during T Q t s . The agreement of both 
estimates of the accidental background of coincidences is well within the statistical 
uncertainties for the detector pairs. In addition, the observed coincidences, n c , 
correspond to both estimates of the accidental coincidence background. This implies 
that no excess coincidences were observed. 

Using EqJ^, we also performed a preliminary analysis of the rates of acciden- 
tal coincidences for three and four-fold configurations of the IGEC observatory. 
The most relevant result here is that the statistical significance of three-fold and 
four-fold time concidences among the current IGEC detectors improves by order of 
magnitudes. 



4. Future Plans and Conclusions 

This first IGEC joint analysis has shown that, with the current detector perfor- 
mances and the selected coincidence time window, at least three detectors simul- 
taneously operating can perform an autonomous search for gw bursts with a very 
low false alarm rate even at signal-to-noise ratio as low as 3 — 5 in single detectors. 
Therefore, the continuation of the IGEC international effort is very strongly moti- 
vated. Moreover, the joint work of the participating groups has instigated efforts 
to co-ordinate data analysis techniques between the different groups. 

No coincidence above the expected accidental background were found in this 
preliminary analysis. Work is in progress to complete the analysis to three and four- 
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fold coincidences, as well as to set upper limits on the rate of incoming gw bursts 
and on the amplitude of single gw bursts associated with selected time windows of 
astrophysical interest. 

Improvements in the detector performances will lead to increase the sensitivity 
of the IGEC observatory in two respects. On one side the thresholds for gw burst 
search will be lowered without increasing the level of the accidental background rate. 
On the other side the effective bandwidths of the detectors will be also widened, 
thus decreasing the uncertainty in the estimated arrival time of a gw burst. The 
latter will allow the lowering of the rate of accidental coincidences and, above all, 
it will give the opportunity to measure the propagation speed and direction of the 
incoming gw. 

In this framework, we think that the participation within the IGEC observatory 
of the interferometric detectors as they will begin observations would constitute a 
very important stage towards the establishment of the future worldwide observatory 
for gravitational waves. 
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